Constraining the dark fluid 
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Cosmological observations are normally fit under the assumption that the dark sector can be 
decomposed into dark matter and dark energy components. However, as long as the probes remain 
purely gravitational, there is no unique decomposition and observations can only constrain a single 
dark fluid; this is known as the dark degeneracy. We use observations to directly constrain this 
dark fluid in a model-independent way, demonstrating in particular that the data cannot be fit by 
a dark fluid with a single constant equation of state. Parameterizing the dark fluid equation of 
state by a variety of polynomials in the scale factor a, we use current kinematical data to constrain 
the parameters. While the simplest interpretation of the dark fluid remains that it is comprised of 
separate dark matter and cosmological constant contributions, our results cover other model types 
including unified dark energy /matter scenarios. 

PACS numbers: 95.36.x,98.80.-k 



I. INTRODUCTION 

The standard cosmological model appeals to two sep- 
arate dark components — dark matter and dark energy 
— and the usual application of observational constraints 
places limits on each of these. However, at present these 
two components have only been detected through their 
gravitational influence, and these measurements do not 
provide enough information to permit a unique decom- 
position into these components. Rather, it is a model 
assumption that the two components are separate. This 
point was first made in Ref. pL,] and specific consequences 
explored in Refs. 0, [!]■ Ref- 0] extended it to coupled 
models and named it the dark degeneracy. 

This degeneracy is extremely general. One can write 
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where G^i, is the Einstein tensor, T^i, the energy- 
momentum tensor, and k the gravitational coupling. Any 
gravitational probe of the dark sector amounts to evalu- 
ating the left-hand side of this equation, and interpret- 
ting the non-zero result as evidence for the dark sector. 
Having obtained the dark sector energy-momentum ten- 
sor this way, such observations can offer no guidance on 
how that tensor might be split amongst different dark 
components.^ 

Although the degeneracy holds even for structure for- 
mation probes (both linear and non-linear) of the dark 
sector, we will focus here only on kinematical probes, i.e. 
those referring to the homogeneous and isotropic back- 
ground cosmology. A complete dark sector description is 



^ We work entirely in Einsteinian gravity. A somewhat related is- 
sue arises in modified gravity models, where the 'non-Einsteinian' 
gravitational terms could be shifted to the right-hand side of the 
generalized Einstein equation and potentially reinterpretted as 
matter terms, see e.g. Ref. [^. 



then given by the total dark sector energy density and its 
equation of state Wdark, the latter determining the evo- 
lution of the former. This enables a simple analysis. To 
include structure formation data, while retaining statisti- 
cal homogeneity, one would also need to consider at least 
the dark sector sound speed and perhaps also anisotropic 
stress [ll, |3| , and one might expect that structure forma- 
tion observations would end up mostly constraining their 
form rather than imposing further upon Wdark- 

As a simple example [Ij, in the standard cosmologi- 
cal model the redshift evolution of the total dark sector 
equation of state 



Wdark 
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(1) (where 'i' runs over the dark components) is given by 
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l-H0.31(H-z)3 



(3) 
(4) 



Here and Qh are the total matter density parame- 
ter and the baryon density parameter, and the subscript 
indicates present value. The second line follows from 
inserting the values rJm.o — 0.274 and Ob,o — 0.046 ob- 
tained from current data compilations Q . 

Inclusion of a single dark energy component with this 
equation of state evolution would give the same observa- 
tional predictions as the standard cosmology. Indeed, 
once one allows the dark energy equation of state to 
evolve arbitrarily, one cannot say anything from obser- 
vations about the dark matter density Qdm, as its effects 
can always be reinterpretted as due to the dark energy. 
Analyses which appear to measure Qdm accurately only 
manage to do so because the dark energy parameteriza- 
tion adopted is not general enough to be able to mimic 
the form of Eq. Q. 
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This degeneracy is perfect. All we can say in a model- 
independent way is that the present total dark sector 
density is about 0.95, and that its equation of state is 
constrained to evolve in a particular way from an early- 
time value at or near zero to arrive at its present value 
of w'^f.^.{z = 0) ~ —0.8. Our aim in this paper is to 
more precisely quantify these constraints, by considering 
general dark sector models that do not include a pure 
cold dark matter component. 



II. MODELS AND DATA 

Having set up this non-standard framework for analyz- 
ing the dark sector, our analysis procedure is standard 
and straightforward. 



A. Models 

In order to impose constraints on the dark fluid, we 
need to employ a parameterization of its equation of 
state. Henceforth we drop the subscript 'dark', using 
w throughout as the total dark sector equation of state. 
As we wish to reach high redshift, we parameterize w 
as a function of scale factor a, which has the bounded 
domain < a < 1. We simply expand this as a power 
series in a at the present, which we find to be sufficient. 
At linear order this is the well-known CPL parameteriza- 
tion w = wq + {1 — a)wi 0, normally applied to the dark 
energy alone but here referring to the combined dark sec- 
tor. To general order this expression was given, again for 
the dark energy alone, in Ref. |8|, which in actual calcu- 
lation considered the linear and quadratic versions. As 
we are demanding that our parameterization describe the 
entire dark sector, we will explore up to cubic order. Our 
models are hence 

N 

u;(a) = ^w;„(l-a)" (5) 

for different choices of N. 

As we will see, the data strongly demand that 'w{a) is 
close to zero as a ^ 0, i.e. the dark sector is required 
to behave as dark matter at early times (note that this 
is true even without inclusion of any structure forma- 
tion data). Accordingly we will also consider the same 
expansions supplemented with the additional constraint 
w{a — 0) ~ which fixes their highest-order term as 
a function of the others, and with the combined con- 
straints w{a = 0) = and dw/da\^^Q = 0. We will call 
these the 'constrained' and 'doubly-constrained' expan- 
sions respectively. The former can be used from linear 
order upwards, and the latter from quadratic upwards. 

In all cases we choose priors on the expansion coeffi- 
cients which are wide enough that the outcome is deter- 
mined entirely by the data and not the prior. 



The standard rulers measured by the cosmic mi- 
crowave background (CMB) and baryon acoustic oscil- 
lations (BAO) are all fixed by the sound horizon at de- 
coupling (which we fix to be Zdec = 1089), which is de- 
pendent on the baryon density. The standard rulers are 
all distance ratios independent of Hq , and the supernovae 
make no measurement of the Hubble parameter Hq since 
it is degenerate with the unknown absolute normalization 
over which we marginalize. As our data compilation does 
not constrain Hq, it is not able to give an accurate con- 
straint on the total dark sector energy density. However 
the physical baryon density fliyh^ is accurately measured 
to the usual value, and so if supplemented with a direct 
determination of Hq, we would find the expected total 
dark sector energy density rJdark — 0.95. In our analysis 
we allow Q\jh^ and Hq to vary, but constrain them us- 
ing other datasets, and marginalize over them, in effect 
marginalizing over the total dark sector energy density. 

In light of the above, when we quote model parame- 
ter counts they are of the dark sector equation of state, 
and do not include the total dark sector energy density. 
For the general dark sector models, we quote the number 
of parameters specifying w{a). For ACDM, the equiva- 
lent parameter is the relative amount of dark matter and 
dark energy at present (equivalently, the coefficient of the 
(1 -I- z)3 term in Eq. (gl). 



B. Data 

We use a fairly typical combination of kinematical data 
to constrain our models. Standard candle data comes 
from supernova type la luminosity distances, for which 
we use the cut Union supernova sample (with sys- 
tematic errors included), and standard ruler data comes 
from the angular positions of the CMB [ll[ and BAO 
peaks [T^. Note that Ref. fl)] give constraints on the 
scaled distance to recombination R and the angular scale 
of the sound horizon la- These are defined to be 

i?^y^r(zcMB), L^^^^p^. (6) 

Since R is scaled by the physical matter density, and 
so makes assumptions about the separability of the dark 
matter and dark energy, we ignore it in this work. We use 
only the constraints on la, as well as those on fly^h^ and 
the correlations between the two. This still assumes that 
the sound speed of the dark component at high redshift 
is small, in order to avoid an early ISW effect shifting the 
peak position. We also include the SHOES (Supernovae 
and Hq for the Equation of State) measurement of the 
Hubble parameter today. Hq = 74.2 ± 3.6 kms~^Mpc~^ 

i. 
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TABLE I: Parameters (of the dark sector equation of state) 
and best-fit clii-squared for our various models. Tfie con- 
strained models force w{a = 0) = 0, and the doubly- 
constrained ones additionally dw/da\^^Q — 0. 



Model 


Dark sector 
parameters 


Xmin 


ACDM 


1 


31f.9 


Constant w 


1 


391.3 


Linear (CPL) 


2 


312.1 


Constrained linear 


1 


320.5 


Quadratic 


3 


309.8 


Constrained quadratic 


2 


311.9 


Doubly-constrained quadratic 


1 


313.5 


Cubic 


4 


309.6 


Constrained cubic 


3 


311.1 


Doubly-constrained cubic 


2 


311.5 



III. RESULTS 

TableUshows the number of adjustable parameters and 
best-fit clii-squared for most models, including ACDM for 
comparison. The total number of data points is 313 (308 
SN-Ia, 2 BAO, 2 CMB and 1 from the SHOES project), 
but correlations between the data points make it difficult 
to state the number of independent data points, and so 
the number of degrees of freedom. We can say that the 
number of degrees of freedom is, at most, 311 minus the 
number of dark sector parameters, meaning that the data 
is an acceptable fit to all the models, except the constant 
w model, if the correlations are small. 

The corresponding w{a) curves of the best-fitting ver- 
sion of each model are shown in Fig. [TJ The immediate 
conclusions from these are as follows: 

1. ACDM, as expected, does a good job of fitting the 
data, bettered only by other models with more dark 
sector parameters. For our data the best-fit flm is 
0.26 ± 0.02. This agrees well with the result of Ko- 
matsu et al. [6], with somewhat larger uncertainty 
as less data is being used. 

2. Even the best-fit version of the constant w model 
is a very poor fit. The dark sector has not had a 
constant equation of state throughout its evolution. 

3. The full four-parameter cubic does not significantly 
improve the fit over the three-parameter quadratic, 
indicating that three-parameter models saturate 
the constraining power of the data. 

4. All the models have w{a = 0) at or very close to 
zero, enforced entirely by the data. This indicates 
that the full phenomenology can be captured using 
the constrained versions of the expansion, reducing 
the variable parameter set by at least one. 
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FIG. 1: The best-fit w{a) for our various models. Note that 
the approach to = at a = is determined entirely by the 
data in the unconstrained cases, while being enforced in the 
constrained models. 




-0.7 -0.65 -0.6 -0.55 




-0.7 -0.65 -0.6 -0.55 0.6 0.65 0.7 

^0 ^1 



FIG. 2: Parameter constraints on the linear (CPL) parame- 
terization. The constrained models lie on the line wi — —wq, 
shown as the (blue) line on the liio, wi plot. 

the data, though not quite as good as the ACDM 
fit. By contrast, the single-parameter constrained 
linear fit is a poor fit to the data when compared to 
the ACDM fit which has the same number of pa- 
rameters (though the general CPL model can ac- 
ceptably fit it by overshooting to w > at early 
times). 

6. In terms of giving good fits to the data for econom- 
ical numbers of parameters, apart from ACDM the 
two-parameter constrained quadratic and doubly- 
constrained cubic fits are the most appealing mod- 
els. 



5. 



The doubly-constrained quadratic (i.e. simply 
w{a) = woo^) gives a tolerable one-parameter fit to 



More relevant than the best-fit parameters are the 
ranges of parameter values permitted by the data. As 
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FIG. 3: Parameter constraints on the quadratic parameteri- 
zation. Here Wb ~ wq + wi + W2, which is the combination 
which is held at zero in the constrained quadratic case. 
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FIG. 4: Parameter constraints on the constrained quadratic 
parameterization, where ■w{a = 0) is set to zero. The doubly- 
constrained models lie on the line wi = —2wo, shown as the 
(blue) line on the wo, wi plot. 



all constant w models are ruled out, we show in Figs. [2] 
and [3] the constraints on the parameters of the CPL and 
quadratic parameterizations, and in Fig. 5] we show the 
constraints on the constrained quadratic model, which 
forces w{a = 0) = 0. 

For illustrative purposes, Fig. [5] shows 400 quadratic 
w{a) curves drawn randomly from the Markov chain. 
They are shaded so that the likelihood increases from 



blue (darker) to red (lighter). We can see that the high- 
redshift region near a = is strongly constrained and 
requires w ~ 0, while the value of the total equation 
of state parameter today, wq = w{a = 1), is only very 
weakly constrained by the data sets used in this work.^ 
Due to the integrated nature of the distance constraints 
on w, curves that oscillate around the best-fit incur only 
a small penalty. This becomes more problematic when 
going to higher order in power law expansions since then 
more oscillations become possible. As the oscillations 
have to average out, the expansion parameters are highly 
correlated and not really independent, which can already 
be seen for wq and wi in the quadratic case in Fig. [31 

We also carried out a different analysis where the equa- 
tion of state w is allowed to take different values in binned 
regions of scale factor a, for simplicity taken to be lin- 
early spaced with 10 or 20 bins. A Principal Component 
Analysis showed that three modes were well measured 
(cr < 0.1), supporting the conclusions above. This analy- 
sis also showed that the best measured mode was peaked 
in the highest-redshift bin. 



IV. CONCLUSIONS 

Due to the dark degeneracy, there is no unique split 
into dark matter and dark energy. For this reason, we 
considered in this paper the total dark sector equation of 
state. We parameterized it with polynomial expansions 
in the scale factor a and used type la supernovae, baryon 



^ As we were completing this work, Mortonson et al. Il3ll arXived 
a paper considering this specific point in much more detail. This 
point had also previously been noted in Ref. \Ti . 
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acoustic oscillations, and the CMB peak position to find 
constraints on the expansion parameters. 

The strongest constraints come at very high redshift, 
from the CMB and BAO measurements of the sound hori- 
zon at decoupling. There the dark fluid must evolve with 
an equation of state close to zero, to recover the correct 
angular scale of the acoustic oscillations. At low redshift 
the constraints from SN-Ia and BAO are weaker, but re- 
quire a negative pressure fluid with Wdark — —0.8. 

We found that ACDM gave the best one-parameter 
fit to current data, bettered only by other models with 
at least two parameters such as the doubly-constrained 
cubic expansion. However, none did very much better 
than ACDM. Accordingly, ACDM remains the most com- 
pelling interpretation of present data, despite the dark 
degeneracy. Nevertheless our polynomial fits indicate the 
constraints that can be applied to more general types of 
dark sector model. 

One could speculate about how future measurements of 
the same general type might tighten our constraints, but 
as remarked in the introduction they cannot address the 
issue of separability of the dark sector into components. 
Rather, to defeat the dark degeneracy what is needed 



are non- gravitational probes of the sector, for instance 
direct detection of dark matterparticles at accelerators or 
in underground experiments [J|. Such detections would 
immediately move such particles from the dark sector to 
the 'known' sector, removing them from the dark sector 
degeneracy. 
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